Soft two-meson-exchange nucleon-nucleon potentials. 
II. One-pair and two-pair diagrams 



Th.A. Rijken 

Institute for Theoretical Physics, University of Nijmegen, Nijmegen, The Netherlands 



V.G.J. Stoks* 

Department of Physics, The Flinders University of South Australia, Bedford Park, South 
On , Australia 5042, Australia 

(version of: February 9, 2008) 



in 



0^ 



X 



Abstract 



■ Two-meson-exchange nucleon-nucleon potentials are derived where either one 
^ . or both nucleons contains a pair vertex. Physically, the meson-pair vertices 

, are meant to describe in an effective way (part of) the effects of heavy-meson 

I exchange and meson-nucleon resonances. From the point of view of "duality," 

' these two kinds of contribution are roughly equivalent. The various possibili- 

■ ties for meson pairs coupling to the nucleon are inspired by the chiral-invariant 



phenomenological Lagrangians that have appeared in the literature. The cou- 
pling constants are fixed using the linear a model. We show that the inclusion 
of these two-meson exchanges gives a significant improvement over a potential 



^ ' model including only the standard one-boson exchanges. 
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I. INTRODUCTION 



In this second paper on two-meson-exchange potentials, we derive the contributions to 
the nucleon-nucleon potentials when either one or both nucleons contains a pair vertex. The 
corresponding "seagull" diagrams are labeled as one-pair and two-pair diagrams, and are 
depicted in Fig. |I[ They are in a different class than the planar and crossed-box diagrams. 
The latter can be understood as the second-order contributions in a series expansion of 
multi-meson exchanges and were calculated in the previous paper |^. The two types of 
two-meson-exchange potentials presented here and in the previous paper |l|] are part of our 
program to study the effect of two-meson exchanges in potential models, and to extend the 
Nijmegen soft-core one-boson-exchange potential p| to arrive at a new extended soft-core 
nucleon-nucleon model, hereafter referred to as the ESC potential. 

Our motivation for deriving the pair- meson potentials comes from "duality" P,^. Since 
we do not explicitly include any contributions to the potentials coming from intermediate 
states with nucleon resonances, we have indicated in Fig. ^ how these resonances can be 
included via meson-pair contributions. According to "duality," the resonance contributions 
to the various meson-nucleon amplitudes can be described approximately by heavy-meson 
exchanges. Treating the heavy-meson propagators as constants, which should be adequate 
at low energies, leads directly to pair-meson exchanges. Hence, the pair-meson potentials 
can be viewed as the result of integrating out the heavy-meson and resonance degrees of 
freedom. We are less radical than Weinberg (see, e.g., Ref. |^), in that we do not integrate 
out the degrees of freedom of the mesons with masses below 1 GeV. 

The pairs we consider are tttt, 7rp, ne, and ee. They are inspired by the chiral-invariant 
phenomenological Lagrangians that have appeared in the literature, see for example 
Here the £(760) scalar meson is treated as a very broad meson with a width of F^ ~ 640 
MeV p] . The potential due to the exchange of a broad meson can be approximated by the 
sum of two potentials where each potential is due to the exchange of a stable meson P,p!0 



Because of the large width of the e(760), the low-mass pole in the two-pole approximation 
is rather small (~ 500 MeV). Hence, it is expected to contribute significantly to the two- 
meson potentials, which is the reason why we explicitly included the ee pair potential. To 
emphasize the importance of the low-mass contribution, we will henceforth denote these 
scalar contributions to the pair potentials by ira and aa, where cr stands for the low-mass 
pole. 

We neglect the contributions from the negative-energy states of the nucleons. We assume 
that at low energies a strong "nucleon-pair suppression" mechanism (i.e., suppression of the 
nucleon-antinucleon NN pairs) is operative. This is due to the compositeness and the large 
mass of the nucleons. Obviously, this is only valid in the low-energy regime. From relativity 
we know that nucleon-pair suppression cannot be absolute. However, we mention that it 
is very plausible from the point of view |]TT| of the nonrelativistic quark model (NRQM). 
Also, in the chiral quark-model picture |jl2| it is rather unlikely that NN pairs play a role 
in the low-energy region. Moreover, nucleon-pair suppression is substantiated by studies 
of the large Nc limit in QCD [|1^. In this case then, the Thomson limit is provided by 
the negative-energy contributions of the constituents, in particular the current quarks. A 
covariant description of this nucleon-pair suppression could very well be done by introducing 
off-mass-shell factors to the meson-nucleon-nucleon vertices. Instead of attempting such a 
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covariant description in this paper, we simply neglect the transitions to the negative-energy 
states. 

The paper is organized as follows. In Sec. |I| the meson-pair-exchange kernels are de- 
rived. We give the interaction Hamiltonians and the vertices in Pauli-spinor space, while 
the implementation of the Gaussian form factors is briefly indicated. In Sec. |IT1| we derive 
the one-pair and two-pair potentials. We then extend the calculations by including the 1/M 
corrections from the pseudovector vertex and from the nonadiabatic expansion of the energy 
denominators. The strength of the pair interactions can be estimated using the saturation 
with one-boson exchanges, but this involves some ambiguity in how to approximate the 
intermediate-boson propagator. In Sec. [V| we therefore choose to fix the pair coupling con- 
stants at their values as obtained in the linear a model [Q. These values do not have to be 



taken exactly, of course, but they clearly suffice to demonstrate the improvements that can 
be obtained when we include these pair interactions in a nucleon-nucleon potential model. 
Plots of the various pair potentials are then shown and discussed in Sec. |VI[ In this section 
we also compare this new ESC potential with the one-boson-exchange Nijm93 potential [jl5 



demonstrating the significant improvement in the description of the scattering data that has 
been obtained by including the two-meson exchanges. 

Finally, some details on parts of the calculation are collected in two appendices, while 
the coordinate-space pair potentials are given explicitly in Appendix |C[ 

II. THE PAIR-MESON-EXCHANGE KERNEL 

The fourth-order two-meson-exchange kernel is derived following the procedure as dis- 
cussed in the previous paper [0, which closely follows an earlier publication [|1^], where we 



derived the soft two-pion-exchange potential. For details and definitions we refer to these 
two references. 

Rather than repeating the derivation again, we here immediately proceed to give the 
one-pair and two-pair kernels. Introducing subscripts IP and 2P, the one-pair kernel is 
given by 

ifip(p', v\W)a'y;al> = -{2Ti)-^[W - W(p')] [W - W(p)] 

^ Y d,p'o dpo / dkio / dk2o / dki / rfka z(27r)"^5^(p' - p - ki - k2) 

qIi yi ^ ^ ^ ^ 

X {[F,, ,(p',p'o; Y>.V^W,F^\-V ~ ki, -Po - 

+ [F,F^;^^(p + ki,po + fcio)r.]('^")[F,„(-p', -p, -po)]} 
X 'F\^\^,p,)F%\-^-p,)\\k\-m\^ibY\ (2.1) 
whereas the two-pair kernel is given by 

^2p(p', v\W)a'u,ah = -{27r)-^[W - >V(p')] [W - W(p)] 

X J dpQ J dpo J dkio J dk2o J rfki J dk2 i{27r)~^S^{p' — p — ki — k2) 
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X {[rj,i(p',Po; p,po)][rj,i(-p', -Po; -p, -po)]} 



X 



F^w\v.Po)F^w\-V.-Po) '[e,-mi + ^5]-\ (2.2) 



-1 



Here mi and m2 denote the two meson masses, Fj and Fj denote the nucleon-nucleon-meson 
vertices, and Fj j denotes the nucleon-nucleon-meson-meson vertex; they follow from the 
interaction Hamiltonians (see below). Because we only consider nucleons in the intermediate 
state, we have a = a' = a" = N and b = b' = b" = N. Note that the first term between 
the curly brackets in the one-pair kernel corresponds to having the pair vertex on the first 
nucleon (hence, no label a"), and the second term to having the pair vertex on the other 
nucleon (no label b"). 

From the explicit expressions ( p.l|) and (|2.2|) , it is clear that one can perform the integra- 
tion over the energy variables p'q, po, fcio, and /c2o- The execution of these integrals is quite 
similar to those worked-out explicitly in Ref. 0], and so are the results. Details are given 
in Appendix In the case of the one-pair diagrams, we have taken the diagram where the 
meson-pair vertex is on line a. 

In the adiabatic approximation, i.e., E{p) ^ M, the energy denominators of the various 
time-ordered diagrams are 

1 1 



1 



1 



2UiUJ2 UJi{lUi + UJ2) 



D^'\u;,,u;2) = -^ (2.3) 

Here we have labeled the three time-ordered one-pair graphs by a, b, c, which correspond, 
respectively, to the diagrams depicted in Fig. |^. However, since we always have to add the 
contribution from the "mirror" graphs anyway (i.e., with the pair vertex at line b), we have 
already included the resulting factor of 2 in the D^^^ energy denominators of Eq. ( |2.3| ). The 
two time-ordered graphs for the two-pair diagrams are each other's "mirror" graphs. In 
most cases the vertices do not explicitly depend on the time ordering and we can add the 
three time-ordered one-pair energy denominators to give 

D«(c.i,c.2) = /^W+/^W + DW = 4^. (2.4) 

UJ1UJ2 

Before we can proceed and calculate the pair-meson potentials, we have to define the 
nucleon-nucleon-meson (NNm) and nucleon-nucleon-meson-meson {NNmim2) Hamiltoni- 
ans. For point couplings the nucleon-nucleon-meson Hamiltonians are 

Up = I^^jj^,^ r^p-d^cf^p, (2.5a) 
T-ts = gsi>r^p-(j)s , (2.5c) 
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where (f) denotes the pseudovector-, vector-, and scalar-meson field, respectively. For the 
isospin 7 = mesons, the isospin Pauli matrices, r, are absent. 

For the phenomenological meson-pair interactions, the Hamiltonians are 



Us = ipip [5((^^)„7r-7r + g(aa)cr'^\ /m^, 



/(tttt)! 

' 2M 



(2.6a) 

(2.6b) 

(2.6c) 
(2.6d) 



The transition from Dirac spinors to Pauli spinors is reviewed in Appendix C of [ITB 



Following this reference and keeping only terms up to order 1/M, we find that the vertex 
operators in Pauli-spinor space for the NNm vertices are given by 







■ fp 




;p) 






cr 










')r«n 


;p) 


= 9v 








» 


= 9s, 







2M' 



2M 



'p' + p) + i(l + Ky)cri X k I ■ 



(2.7a) 

(2.7b) 
(2.7c) 



where we defined k = p' — p and Ky = fv/9v- In the pseudovector vertex, the upper (lower) 
sign stands for creation (absorption) of the pion at the vertex. Similarly, the NNmim2 
vertices result in 



«(p')r?)«(p) 



9{mr)o/n^n Slid 9{(Ta) 



^5'(7r7r)i 

9i7Tp)l 



(±^1 ^^2) + J^[^'0^^ ~ ^^2) - (1 + Ki)cri-(ki X ks 



cri-(ki - k2 



M 



cri-q(±u;i ip UJ2) 



/ml, 



(2.8a) 

/ml, 
(2.8b) 

(2.8c) 
(2.8d) 



where q = |(p' + p) and Hi = (//(7)(^^)j. Again, the upper (lower) sign in front of Ui 
and UJ2 refers to creation (absorption) of the meson at the vertex. For both the NNm and 
NNmim2 vertices, the expressions for «(— p')rM(— p) are trivially obtained by substituting 
(p',p,ki,u;i,cri) (-p', -p, -kj, cj^, erg). 

The generalization of the interaction kernels to the case with a Gaussian (or any other) 
form factor has been treated and explained in ||T6[. We make the substitution 



[k^ -m^ + iS\ 



-1 



k? — jj? + id' 



(2.9) 



for each meson-exchange line in the Feynman diagrams. Here, p(/i^) is the spectral function, 
representing the form factors involved in meson exchange. At low and medium energy, we 
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have to a very good approximation t = k"^ — < 0, and so for space-like momentum 
transfers we can use Gaussian form factors -F(k^) = exp(— k^/A^), where A denotes the 
cutoff mass. The Gaussian form factor is introduced by the substitution 

= ^ _ _F(k5_ 

The NNm and NNmim2 vertices have different form factors. We will use 

i^ivivmim2(ki,k2) = exp(-k2/2A2) exp(-k2/2A2), 

F^ivm(k') = exp(-kV2A^), (2.11) 

where Ai and A2 are the form factor masses for mesons mi and m2, respectively. A moti- 
vation for this prescription could be that in "duality" the structure of the NNmim2 vertex 
is either saturated by heavy mesons or meson-nucleon resonances. In this last case, assum- 
ing that the meson-nucleon resonance transitions all have roughly the same (inelastic) form 



factor, the form (2.11) is a natural one. 



III. MESON-PAIR POTENTIALS 

In this paper we restrict ourselves mainly to the pion-meson pair potentials. The only 
exception is the cxa-pair potential, but the derivation of the potential for any other combi- 
nation of two mesons will be straight forward. Here a stands for the low-mass pole in the 
two-pole approximation [||,|10| of the broad e(760) scalar meson. (At this point it might be 
worthwhile to mention that in the previous paper |[l[] we refrained from evaluating the aa 
planar and crossed-box contributions. The reason for this is that the leading-order contribu- 
tion of a potential due to the exchange of two isoscalar-scalar mesons is identically zero.) In 
the description of the pion-meson pair potentials we always assign the index 1 to the pion, 
and the index 2 to the other meson, which can be a pion as well. 

Armed with the Pauli-spinor vertices given in the previous section, it is now straight 
forward to derive the one-pair and two-pair potentials. They can be succinctly written as 



p,,,,^^, ^ '■^^^ ^ / / (27r)6 



(3.1) 



where the index n distinguishes one-pair {n = 1) and two-pair (n = 2) meson-pair exchange, 
and (a/3) refers to the particular meson pair that is being exchanged. The product of the 
coupling constants in the two cases is given by 

5,(1) (a/5) = g(^aP)9NNa9NN^, 

9^'\aP) = gl,^, (3.2) 
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with appropriate powers of m^r, depending on the definition of the Hamiltonians. The 
energy denominators D^'^ are given in Eq. (|2.3| ), with the index p labehng the different 

time-ordered processes. Finally, the momentum-dependent operators O^^f^p are given in 
Tables | and [I^. For completeness, these Tables also contain the isospin factors C(")(a/?) as 
derived in Appendix The momentum operators for (7r7r)o and (tttt)! both contain a term 
antisymmetric in ki ^ k2. They only contribute when we make the nonadiabatic expansion 
(see Sec. In the leading-order potential (|3.1|) they drop out when we integrate over ki 
and k2. 

The time-ordering label p in Eq. is only of interest for the (vrvr)! one-pair potential, 
where we have an explicit (101,102) dependence in O^^^^. For all other one-pair potentials we 
have the same energy denominator D^^\uji,uj2) as given in Eq. (|2.4|) , and we can drop the 
p index in O^^^p. For (a/3) = (tttt)!, the three time-ordered diagrams contribute as 

2 



UJlUJ2{uJi + UJ2) 



Similarly, the two-pair potentials are all seen to have the energy denominator D*^^^ of 
Eq. (^.31), where the additional (cui — UJ2Y dependence in O^^j for (vrvr)! can be rewritten as 

{u;,-uj2fD^'^ = ^ ^ ^ 



2uJi 2U02 OOi + 002 

The evaluation of the momentum integrations can now readily be performed using the 
methods given in p],p!6|. There it was shown that the full separation of the ki and k2 depen- 



dence of the Fourier integrals can be achieved in all cases using the A-integral representation. 
Starting-out from Eq. (|0| ), this procedure gives the following generic form of the potentials 
in coordinate space: 

V;S(a/?) = C(")(a/5)(7(")(a/?) Jm ^ 

(3.3) 

The different functions Ba/siri, r2) that occur in this expression involve the functions l2{m, r) 



as defined in Ref. [|T6[, /o(A,r) as defined in Ref. and the integrals 

2 roc 



Bofi{ma,ri;mfs,r2) = - / c/A A^Fa(A, ri)F^(A, r2), 

77 JO 

Bi^i{ma,ri;mf3,r2) = - / dXFa{X,ri)Ffs{X,r2), (3.4) 

TT JO 



F„(A,r) = e-" /^«/2(Vm2 + A^, r). (3.5) 



where 



The differentiation operations needed in Eq. ( |3.3| ) are listed in Appendix A of the pre- 
vious paper [Q, and the resulting coordinate-space potentials are here given explicitly in 
Appendix 0. 
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IV. 1/M CORRECTIONS 

The nonadiabatic correction from the 1/M expansion of the energy denominators is 
explained in Ref. |jT6|. The expansion of the energy denominator involves a momentum 
dependence which can be rewritten in the form [ki-k2 — q-(ki — k2)]/2M. 

Taking out the momentum-dependent factor, the one-pair energy denominators of 
Eq. ( p.3|) give rise to the energy denominators 



2uJiUJ2 t^K^i + ^2) ' 

1/1 1 

+ 



1 1 



(4.1) 



2u;iCo'2 ujiiuji + UJ2) 

Again, the time ordering is only of importance in the (vrvr)! potential, where the contributions 
sum up to 

2 



r-jna 



(4.2) 



In all other cases the three energy denominators of Eq. ( [4.1D can be summed directly, yielding 



1^1 



Ul UJ2 UJ1+UJ2. 

which gives rise to a coordinate-space function of the form 



(4.3) 



BTM,T2) 



2 r°°d\ 



TT Jo 



I2{ma,ri)l2{mp,r2 
-Fa{X,n)Fp{X,r2) 



(4.4) 



The energy denominators D'^^{uji,uj2) are symmetric under interchanges of labels 1 and 
2, and so the term proportional to q-(ki — k2) will only contribute in combination with the 
antisymmetric terms in 0^^\'ki,uJi;'k2UJ2); i.e., in (7r7r)o and (7r7r)i. We find 

2 



V^-[(vrvr)o] = (^] ^ //^^^|^ e^('^^+'^^)-F.(k?)F.(k^) 



krk2)^ + iq-(ki-k2)(o-i + o-2)-(ki x ks 



D°-(^i,^2) 



(4.5) 



[rrT2 

mi 



M 



(2vr) 



X 



^ki-k2)2 + |q-(ki - k2)(<Ti + 0-2)-(ki X ks 



2 2 ' 



y-[(7ra)] 



9{(7cr) dNNo 



M 



(27r)6 



F^{kl)F,{kl){k,-k2)D^^{u;,,u;2) 



(4.6) 
(4.7) 



(Ti-T2) 



x(ki-k2) 



9i7Ta) fNNiT dNNa f f d^kid^h 



(27r) 



!2^.(ki+k2).r^^^j^2)^^^j^: 



cri-kicr2-ki - i(cri ■kicr2-k2 + cri-k2cr2-ki) 



D-^{u;,,u;2), (4.8) 
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where we substituted the isospin dependence and the couphng constants. 

The pseudovector vertex gives rise to 1 /M terms as shown in Eq. ( p. 7a ) . Taking into ac- 
count the different time orderings and discarding the terms antisymmetric under interchange 
of ki k2, we find the following contributions 



vrvr 



9{n7T)o ( InNtt 



\ m. 
,2 , 1,2 



3 

M 



(27r)6 

k^ + k^-z(cri + (T2)-(ki + 



g^(ki+k,).r^^(j^2)^^^]^2) 

k2) X q"' ^ 



I InNtt 



1 

M 



X 



k2 



g{-Ku) InNtt 9NNcr 

ml M 



<fkid^k2 
(27r)6 



'ki , k2\ 

— + — xq 



g^(k,+k,).r^^^]^2)^^(]^2) 



lTrT2 



<fki<fk2 
(27r)6 



X 



cri-k2cr2-k2 - i(cri kicr2 k2 + cri k2cr2 ki 



e'('^^+'^^)"F,(k?)F,(k 
1 



UJlUJ2{uJi + UJ2) 



(4.9) 



(4.10) 



(4.11) 



V. THEORETICAL CONSTRAINTS 
ON PAIR COUPLING CONSTANTS 

In principle, each of the pair-meson potentials contains a free parameter, the pair-meson 
coupling constant. A possible way to fix the values for the pair-meson coupling constants 
is to assume that the coupling of a meson pair {a(3) to a nucleon is dominated by an 
intermediate boson H . We can then relate the pair coupling constant g^ap) to the meson- 
decay and meson- nucleon coupling constants, guap and gNNH, respectively. The relevant 
nucleon-nucleon-meson and meson-decay Lagrangians can then be obtained by extending 



the linear a model with vector and axial-vector mesons; see, for example, Refs. |1T7|-[19 

Meson saturation is graphically presented in Fig. ^. Because of the additional meson 
propagator and the meson-decay vertex, we find the relation 

g{af3) ~ —gNNHgHap, (5.1) 



m 



H 

with appropriate powers of meson masses, depending on the actual form of the meson-pair 
and meson-decay Hamiltonians. Here S denotes the spin of the intermediate meson H and 
we assumed that is much larger than the momentum-transfer squared. However, for the 
(7r7r)i and (vr7r)o pairs the decay can proceed via the p and a mesons, respectively. These 
are substantially lighter than 1 GeV, and so we are faced with the difficult task of finding 
a suitable average momentum squared to approximate (k^ +m%). Another problem is that 
a specific meson pair might be the decay product of a range of different mesons, and so 
Eq. ( p.l[ ) becomes a sum over mesons Hi, confusing the issue even further. Alternatively, 
we can also use the linear a model to generate the pair vertices explicitly, which we shall do 
in the following. 
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The a model has been discussed extensively in the hterature, and here we only briefly 
outhne its contents to define the quantities we need. The model contains an isotriplet of 
pseudoscalars, tt, and an isosinglet scalar, cr, grouped into 

E = (7-iT-7r, (5.2) 

which transforms under global SU(2)j;,xSU(2)ij as 

E ^ LEi?^ (5.3) 

The nucleon wave function has left and right components, ipL,R = ■, transforming 

as 

il)L LipL, i/jR RipR. (5.4) 
The transformations can be made local by introducing the left and right gauge fields, 

i^^l-r-ip^ + a^), (5.5) 



and their field strength tensors 

li^i^ = di^lv - dj^ + igv[ln, (5.6) 

and similarly for r^j,. The vector and axial vector mesons are given a mass by introducing 
a chiral-symmetry breaking mass term 

£^=imjTr(/^P + rX)- (5-7) 

After the introduction of the gauge fields, the chiral symmetry can be restored by defining 
the following covariant derivatives for the nucleon and (cr, tt) fields, 

^mV' = (^M + havr-Pii + ^gvr-a^j5)i^, 

T^ix^ = du.T, + igvln^ - igv^r^,. (5.8) 
A possible chiral-invariant Lagrangian is now given by 

+g2{i^LT.T.'^T.ijR + ^fiEtEEtV'i) 
+^C(7A^E7^2)^Et^^ + V5^Et7^I?^E^^) 
-iTr(/^,/^'^ + r,,r^^) + iTr(P'^Etp^E) 

-i//2Tr(EtE) - |A2Tr(E^E)2, (5.9) 

where gi, g2, C, fj., and A are free parameters. 

The symmetry can be spontaneously broken by adding a term linear in the a field, 

^SB = Umla, (5.10) 
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with /,r = 92.4 MeV the pion decay constant. Choosing the ground state as (a) = f , this 
spontaneous symmetry breaking introduces a mixing between the tt and fields, which 
requires a redefinition of the fields, given by 

cr ^ cr + W, 

TT^^TT, (5.11) 

rup 

where we defined m\ = rn?^ + {gyvY, v = (mA/mp)/.,,, and D^^tt = d^iz + gyi^ x p^. The 
shift of the a field gives the nucleons a mass M = {gi — g2v'^)v, which imposes a constraint 
on 5(1 and g2. 

Comparing the various pieces of the full Lagrangian, Co + Cm + ^^sb, with the interaction 
Hamiltonians of Sec. |T|, we identify the following relations for the single-meson coupling 
constants: 

- 3g2V^) = M/v - 2g2v'', 




9nna — /awtt- (5.12) 



The numerical value for C is obtained by imposing the Goldberger-Treiman relation, 

/aWtt gA /r; 1Q\ 

2/^ 

with gA = 1.2573 the weak interaction axial-vector coupling constant. Similarly, the meson- 
pair coupling constants are found to be 

/ M ''tia\ 

9{TnT)o = ~ TTT 7 9NN(7 



9{aa) 



9(7T7t) 




nip 



9NNp, , -.x 

9(^Ph = Y^{9A + l), 

2/2 niAmp rn\J ' 

Working-out the various pieces of the Lagrangian we can also find relations for the 
meson-decay coupling constants. It is then found that the relation ( |5.1| ) is indeed satisfied, 
although in some cases the approximation is rather crude, as we already suspected. However, 
the meson-saturation assumption is still very useful to suggest the obvious relation 

/(tttt)! = 9(Tnr)i {fNNp/ 9NNp)- (5.15) 
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VI. RESULTS AND DISCUSSION 



The complete one-pair and two-pair potential can be written as 

V{a[3) = V^^\a(3) + V^^\a[3) 

+\/""(a/3) + \/P"(a/3), (6.1) 

for the pairs (a/3) discussed in this paper. To present the potentials in graphical form, we 
employ the meson-nucleon coupling constants and cutoff masses of a preliminary version of 
the Nijmegen extended soft-core (ESC) potential. The ESC potential is still under construc- 
tion, but the present values for the coupling constants clearly suffice for illustrative purposes. 
The meson-nucleon coupling constants and cutoff masses are already listed in Table IV of 
the previous paper , while the pair coupling constants are here given in Table |T|. All pair 
coupling constants are fixed at their theoretical values as given in Sec. 0. This constraint 
might be a bit too severe, but the values can at least be expected to be reasonable. Further- 
more, in this way we can get a feeling about the importance of two-meson exchanges with 
respect to one-boson exchanges, because we do not introduce any new free parameters. 

In Figs. ^ and we compare the various types of one-pair and two-pair exchanges for 
/ = and / = 1, respectively. They consist of scalar 0"'""'" exchanges [(7r7r)o and (o"(t)], 
vector 1 exchange [(vrvr)!], and axial l"*"^ exchanges [(vrp)i and (ttct)]. The central 1 = 
potential is clearly dominated by the one-pair (7r7r)i potential, whereas the scalar and vector 
two-pair potentials largely cancel each other. For 7 = 1, the scalar one-pair and two-pair 
potentials are the same, whereas the vector potentials change by a factor -3, resulting in a 
much smaller (but still attractive) central potential. For the spin-spin and tensor potentials 
the 7 = and 7=1 cases differ by a simple factor of -3, where the two-pair axial tensor 
component turns out to be completely negligible. 

In all cases, the 1/M corrections from the nonadiabatic expansion of the energy denom- 
inators and from the pseudovector-vertex correction in the pion vertex function have the 
opposite sign. In the (vr7r)o potential the cancelation is even exact. The 1/M corrections to 
the spin-spin and tensor potentials are only due to the (ttct) potential, whereas the (vr7r)o 
and (7r7r)i potentials only contribute to the central and spin-orbit potentials. 

In Fig. 1^ we compare the spin-orbit components from the one-pair potentials. They are 
all due to the 1/M nonadiabatic or pseudovector-vertex corrections in the (vr7r)o and (7r7r)i 
potentials. The (tttt)! potential also has a spin-orbit potential due to the 1/M term in the 
(7r7r)i pair vertex ( |2.8b| ), as given in Eq. ( P5| ) of Appendix |C[ This contribution is exactly 
the same as the nonadiabatic contribution, which means that in the (vrvr)! potential there is 
an even larger cancelation between the nonadiabatic (long-dashed line) and pseudovector- 
vertex (dash-dotted fine) spin-orbit parts than Fig. ^ might suggest. Note that in the 
(7r7r)o potential the cancelation between the nonadiabatic (solid line) and pseudovector- 
vertex (short-dashed line) spin-orbit parts is exact, and so there is no spin-orbit potential 
from O"*""*" pair exchange. Obviously, if we were to take the pseudoscalar coupling, ipij^ipcf), 
for the pion, rather than the pseudovector coupling, il'lblii'^^d^'Pi these cancelations do not 
occur, and we are left with the spin-orbit contribution from the nonadiabatic expansion. 

Finally, we should mention that in some cases the two-pair diagrams are in principle 
included in the cases where we use the exchange of broad mesons. This is the case for the 
two-pair (7r7r)o and (7771)1 potentials. Here we use in the Nijmegen work a broad e and 
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p one-boson-exchange potential. If such a broad meson exchange is included exactly, the 
two-pair vrvr potentials should be omitted. When such broad mesons are included in a two- 
pole approximation [p|,p!0|, however, a two-pair contribution might still be useful. Whether 
these two-pair contributions in that case should be included in full or partially suppressed 
is presently under investigation. 

Let us now return to a comparison of this new ESC potential with the one-boson-exchange 
Nijm93 potential |TB|. It is important to realize that both models contain essentially the 



same set of parameters, and so the ESC model is a true extension of the Nijm93 model. 
Using the linear a model as a means to get a reasonable estimate for the various meson- 
pair coupling constants, however, the contributions from the two-meson exchange diagrams 
(the planar and crossed-box diagrams as evaluated in the previous paper as well as the 
one-pair and two-pair diagrams as evaluated in the present paper) are included without the 
introduction of any new parameters. The 14 free parameters of the ESC model are fitted to 
the 1993 Nijmegen representation of the hypersurface of the NN scattering data below 



^lab = 350 MeV |2T[, updated with the inclusion of new data which have been published 
since then. 

The results for the 10 energy bins are given in Table where we compare the results 
from the updated partial-wave analysis with the Nijm93 and ESC potentials. Clearly, the 
inclusion of the two-meson exchanges in the ESC model allows for a substantially better 
description of the AW scattering data than what could be achieved with the one-boson- 
exchange Nijm93 model. At present, the ESC potential is the only meson-theoretical model 
which can give such a good description of the scattering data (using only a limited set of free 
parameters). The quality of the ESC potential is even better if we restrict its application to 
energies below Tiab ~ 300 MeV, giving x^/-^data = 1-137 for the 0-290 MeV energy interval. 
This was to be expected, since the nonadiabatic expansion in principle is only valid below 
the pion-production threshold. 

To summarize, the combined results of the present and the previous paper show that it 
is possible to construct a nucleon-nucleon potential model, based on a chiral-symmetric La- 
grangian, that gives a good description of the nucleon-nucleon scattering data. Fine-tuning of 
the 14 parameters of this extended model allows for a substantially better description of the 



data than what was possible with the 14-parameter one-boson-exchange Nijm93 model [15 
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APPENDIX A: ENERGY INTEGRALS 

We discuss here the treatment of the energy integrals that occur in the evaluation of 
the meson-pair potentials. To prevent displaying all the indices, we introduce the following 
convenient notations 
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2, uj' = Vk'2 + m'2, 



A = E(p) - A' = E{p') - \W, A" = E{p 



2 '^'^ ' 



(Al) 



where p' — p = k + k' and p" = p + k. Note that with this notation, k and k' correspond 
to ki and k2, respectively, introduced in the main text, 
(i) The one-pair graph. Here we encounter the integral 

Jip(p', p\W) = -(27r)-2 [W - 2E{p')] [W - 2E{p)] 

X J dpQ J dpQ J dk^ J dko 5^{p' — p — k — k') 



x[k"-m" + ^6]~' El^\p',p',)-' Fl^\-p',-p', 



w y 



p - k, -po - ko 



-1 



x[e-m' + z6]'' El^\p,po) E^\-p,-po 



(A2) 



Similar integrals were treated in |T^], to which we refer for the method of evaluation. We 
first bring the integral into the form 

/+00 r+oo r+oo r+oo r+oo r+oo r+oo 

da dp dp'o / dpo / dk'Q / dko / dp^ 
-oo J—oo J—oo J—oo J—oo J—oo J —oo 

X e*"(Po-'=o-Po')g«/3(Po'-Po-fco) 

X [uj'^ - k'^ - i5]-^[uj^ -kl- i5Y^ 

X [A'^ - p'^ - i5]-^[A" + pI - i5]-^[A^ -pI- i5]-\ (A3) 

The energy-variable integrations can be performed in a straight forward manner using the 
residue theorem, e.g. 



dkn 



c<j2 — fco — i5 2uj 



(A4) 



where in the exponential the (-)-sign and the (+)-sign apply to /? > and P < 0, respectively. 
Also, 



+ 00 



„ e'fPo ^ j 2m e-'f^" , /3 > 
A" + p^' - z5 ~ 1 , /? < ■ 



Keeping track of the signs in the exponentials, the intermediate result is 



(A5) 



oo roc 



da I dp e-'°("'+^'-^")e-*'^("+^+^") 



roo 

da / dp 

-oo Jo 



+ia{Lu'+A'+A") -if3{LU+A+A") 



oo J a 



+ I da I dp ^+i^(^'+A'+A")^+if5(u.+A-A") 



(A6) 



Performing the remaining elementary integrals, we end up with 



14 



Jlp 



27ii 



1 



1 



AujLj' [A' + A" + Lj' A + A' + UJ + Lj' 

1 1 
+ -r. 77. : 77. + 



1 



A' + A" + uj' A + A" + CU A + A" + UJ A + A' + UJ + uj' 



(A7) 



The terms in the curly brackets correspond to the different time-ordered graphs (a), (b), 
and (c) of Fig. ||. 

(ii) The two-pair graph. Here the integral to be performed is 

J2p{p, P\W) = -(27r)-2 [W - 2E{p')] [W - 2E{p)] 

X J dpQ J dpo J dkQ J dko 5'^{p' — p — k — k') 



- + i5]-' f}^\p,Po) F^>{-p,-Po 



(A8) 



This integral can be brought into the form 



(A9) 



/+00 /•+00 r+co r+co r+oo 

da / dp'o / dpo / dk'o / dko e^'^Po-Po-ko-K) 

X - k'^ - i6]-^[u;^ ~kl- i5]-^ 
x[A'~p'^-^5]-'[A'^pl-^5]-\ 

The evaluation of the energy-variable integrals gives in this case 
J^p = -{2n)-'\2my [AculuT' [l^ da e--(-+-'+^+^') + da e+-(-+-'+^+^')| . (AlO) 

The two terms correspond to the time-ordered graph (d) of Fig. |^ and its "mirror" graph. 
Carrying out the final a integrations leads to the final expression 



>J2P 



27Ti 



Aujuj' A + A' + UJ + UJ' 



(All) 



APPENDIX B: ISOSPIN FACTORS 

In this Appendix we review the calculation of the isospin factors. The isospin factor 
comes from the evaluation of all possible contractions for the vacuum matrix element of the 
product of the meson fields coming from the contractions of the interaction Hamiltonians 
at the different vertices. For the one-pair diagrams we label the NNmim2 vertex with (1), 
the meson vertex on the other nucleon line with (2), and the pion vertex with (3). The pion 
is propagating with momentum ki and the other meson with momentum \^2■ The isospin 
factors are then found to be 



15 



(7r7r)o : (0|7r-7r(l)T2-7r(2)T2-7r(3)|0) = {SimSin + SinSrm)T2mT2n = 6, 

(aa) : (0|aa(l)a(2)a(3)|0) = (12, 13 + 13, 12),ontractions = 2, 

(7r7r)i : (0|Tr7r x a^7r(l)T2-7r(2)T2-7r(3)|0) = £,,fcrifer2^r2,(0|7r,(l)9''7r,(l)7r„,(2)7r„(3)|0) 

= £ijkTlkT2mT2n 
= 2tT,-T2id^-d^), 
(7rp)i : (0|Tr7r X p(l)T2-p(2)T2-7r(3)|0) = eijkrikT2mr2nSinSjm = -2iTi-T2, 
(tTO-) : (0|ri-7r(l)o-(l)o-(2)T2-7r(3)|0) = Tu6imT2m = T1-T2. 

For the two-pair diagrams, we label the NNmim2 vertices on each nucleon line by (1) and 
(2), respectively. We then find 

(7r7r)o : (0|7r -77(1)77 -77(2) |0) = {dim^im + 5fm^im) = 6, 

{aa) : (0|cra(l)aa(2)|0) = 2, 

(tttt)! : (0|rr77 x 9'^77(l)r2-77 x a^77(2)|0) = £,,fc£„„irifcr2i(0|7r,(l)a^7r,(l)7r„(2)a'^7r„(2)|0) 

= 2Ti-T2 (92^92^-9^92^) 

= TrT2(9r-92^)(9r-92^), 
(7rp)i : (0|rr77 x p(1)t2-77 x p(2)|0) = eijkemni'rikr2Am5jn = 2ti-T2, 
(vra) : (0|Ti-77(l)a(l)T2-77(2)a(2)|0) = rn5,^r2^ = Ti-T2. 

We should point out that for the (7r7r)i vertices the isospin factors are defined to be 2iTi-T2 
and TrT2 for the one-pair and two-pair potentials, respectively, because the 9f dependence, 
i.e., the kj dependence, is already included in the operators O^'^^'p given in Tables | and 0. 



APPENDIX C: COORDINATE-SPACE POTENTIALS 



The explicit formulas for the soft-core meson-pair potentials are given below, 
(i) Scalar vr vr (/ = 0) pairs: 



(1) _ «5'(7r7r)o JnNtt 



K:^r(r) = 6 



(ii) Scalar a ® a pairs: 



mi TT Jo 



TT JO 



mt 



(CI) 
(C2) 

(C3) 
(C4) 



(iii) Vector vr ® vr (J = 1) pairs: 
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3M 



%(-pAK(A,r)f--il[/^,Jr)]^L.S 
Jo M r"^ ' 



'12 



^p2i(r) = - (T1-T2) TdA F.(A,r) 

TT Jo 

(iv) Axial tt (8) p (/ = 1) pairs: 



A! 



e 



2A2F,(A,r) 



>, (C5) 
(C6) 



V;^ailW = -i(Tl-T2) f^^- ^ 



m-j, Mj 



{1 + K,) I'I'M [5i2-2(tri-tT2)] 



+ (^2,. + hp{r) {<T,-a,) + (/^,, - h',^^) hAr) S^X (C7) 



V^lir) = -(Ti-T2)(o-i-a2)^- TciA F.F,(A,r). 

777-^ TT Jo 

(v) Pseudovector tt cr pairs: 

fl'(7r(T) /aWtt 



2 



(r) (o-i-<72) 



+ 



1 



^2,7r - -^2,n ] ^2,a 



T' T' 



2 



rfA< 



TT Jo 
+ 



F': + -F'\F,+ iF': + -F'\F^-2F'^F, 



(C8) 



(C9) 



F^--F'\F,+ {F':--F'\F^-2F'^F'^ 



(A,r)5i2 . (CIO) 
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TABLES 



TABLE L The one-pair isospin factors C^^^(q/3) and momentum operators O^^'j p(ki, wi; k2, u;2)- 
The index p labels the three time-ordered contributions of Fig. |3| and is only of relevance for the 
(vrvr)! entry, where they are shown as a row vector. Note that ki = (f /g)(TTiT)i- 

{a(5) CW{ap) 0^^j^^(ki,^i;k2,u;2) 



(vrvr)o 6 
(aa) 2 

(vrvr)! 2iTi-T2 



-ki • k2 + |(<Ti + 0-2) • (ki X k2 



kik2 

i 



- §(<Ti + <T2)-(ki X k2) {UJ1-UJ2; -W1-W2; -W1+W2) 

(1 + Ki)cri-(ki X k2)<T2-(ki x k2) + §(cri + cr2)-(ki x k2)q-(ki - k2] 



2iT-\ -To — 

M 



<Ti-kicr2-ki + i(l + Kp)(cri-kicr2-k2 + <Ti-k2cr2-ki - 2<ti •cr2ki •k2^ 



(vrcr) T1-T2 



[cri ki<T2 k2 + cri k2<T2 ki - 2cri kicr2 ki] 



TABLE IL The two-pair isospin factors C^'^\af3) and momentum 
operators O^^^ ^(ki, wi; k2, u;2). The index p labels the two time-ordered contributions of Fig. ^ 
but is of no relevance in the final result. 



(a/3) 


C(2)(a/3) 






k2 


^2) 


(7r7r)o 


6 


1 








(aa) 


2 


1 








(vrvr)! 


T"rr2 


(cjl - UJ2] 










2TrT2 


(T\(T2 








(vro") 


TI-T2 


cri-(ki - 


k2; 


0-2 


(ki - k2) 



TABLE in. Pair- meson coupling constants employed in the potentials shown in Figs. |5| to ^. 
Coupling constants are at k^ =0. All coupling constants are fixed at their theoretical values as 
derived in Sec. [v|. 







(a/3) 


5/4vr 


//4vr 


0++ 




(vrvr)o 


-0.412 




0++ 




(oa) 


-0.482 




\— 




(vrvr)i 


0.058 


0.216 


1++ 




(7rp)i 


0.598 




1++ 




(vro") 


0.053 
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TABLE IV. s-^d P^r datum (Xp j p ) at the 10 energy bins for the updated partial-wave 
analysis (PWA) and the ESC and Nijm93 potential models. A/'jata lists the number of data within 
each energy bin. The bottom line gives the results for the total 0-350 MeV interval. 



PWA ESC Nijm93 



Bin(MeV) 






A.r> H r> 






x' 




0.0-0.5 


145 


144.45 


0.996 


146.15 


1.008 


185.78 


1.28 


0.5-2 


68 


42.97 


0.632 


47.38 


0.697 


55.06 


0.81 


2-8 


110 


106.28 


0.966 


115.19 


1.047 


128.92 


1.17 


8-17 


296 


276.31 


0.933 


326.31 


1.102 


368.20 


1.24 


17-35 


359 


279.54 


0.829 


332.10 


0.925 


393.92 


1.10 


35-75 


585 


567.18 


0.970 


697.19 


1.192 


1337.28 


2.29 


75-125 


399 


409.58 


1.027 


421.89 


1.057 


480.96 


1.20 


125-183 


760 


820.69 


1.080 


936.53 


1.232 


1443.10 


1.90 


183-290 


1047 


1035.48 


0.989 


1261.50 


1.205 


1995.71 


1.91 


290-350 


992 


997.02 


1.005 


1706.75 


1.721 


2866.39 


2.89 


0-350 


4761 


4697.50 


0.987 


5990.99 


1.258 


9255.32 


1.94 
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FIGURES 



FIG. 1. Feynman one-pair and two-pair diagrams. 

FIG. 2. "Duality" picture of the meson-pair potentials. Ri represents a nucleon resonance, Hj 
the intermediating heavy boson, and Mi and M2 are the mesons being exchanged. 

FIG. 3. Time-ordered (a-c) one-pair and (d) two-pair diagrams. The dotted line with momen- 
tum ki refers to the pion and the dashed line with momentum k2 refers to one of the other (vector, 
scalar, or pseudoscalar) mesons. To these we have to add the "mirror" diagrams, where for the 
one-pair diagrams the pair vertex occurs on the other nucleon line. 

FIG. 4. Meson saturation of a pair vertex. The mesons are denoted by Mi and M2, and the 
intermediating heavy boson by H. 

FIG. 5. Central, spin-spin, and tensor components of the / = one-pair and two-pair potentials 
for (a) r < 1 fm and (b) 1 < r < 2 fm. The pairs considered are 0^^: (vr7r)o and (era); 1 : (vr7r)i; 
and 1++: (vr/9)i and (ttct). 

FIG. 6. Same as Fig. |, but for / = 1. 

FIG. 7. Spin-orbit contributions of the one-pair O"*""*" and 1 potentials for both 7 = and 
/ = 1. A distinction is made between the nonadiabatic and pseudovector- vertex contributions. 
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Figure 1, paper II ( Rijken, Stoks ) 
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Figure 2, paper II ( Rijken, Stoks ) 




Figure 3, paper II ( Rijken, Stoks ) 
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Figure 4, paper II ( Rijken, Stoks ) 
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Figure 5, paper II (Rijken, Stoks) 
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Figure 6, paper II (Rijken, Stoks) 




Figure 6, paper II (continued) 




/ 1=0 

/■ nonadiabatic, 0"^"^ 

/' pseudovector, 0"^"^ 

— - nonadiabatic, 1" 
pseudovector, 1 " 



0.25 0.5 



0.75 1 
r(fm) 




nonadiabatic, 
pseudovector, 0"^"^ 
nonadiabatic, 1" 
pseudovector, 1" 



1 .25 1 .5 



1.75 2 

r(fm) 



500 



> 

O 
05 
> 



-500 




nonadiabatic, 0"^"^ 
pseudovector, 0"^"^ 
nonadiabatic, 1" 
pseudovector, 1 " 



0.25 



0.5 



0.75 



10 



> 

o 

05 
> 



-10 




1=1 



nonadiabatic, 
pseudovector, O"^"^ 
nonadiabatic, 1" 
pseudovector, 1" 



1.25 



1.5 



r(fm) 



1.75 2 
r(fm) 



Figure 7, paper II (Rijken, Stoks) 



